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CLAIMS 

1 . A method of determining a distribution of one or more properties in a medium 
illuminated with radiation from one or more sources, the method comprising the steps of: 

a. receiving the radiation exiting the medium; 

b. deriving one or more optical properties of the medium using one or more 
calibration factors and the radiation received in step (a), wherein the 
calibration factors are variables; and 

c. determining the distribution in the medium using the one or more optical 
properties derived in step (b). 

2. The method of claim 1, wherein at least one of the calibration factors is a source 
coupling factor. 

3. The method of claim 1, wherein at least one of the calibration factors is a detector 
coupling factor. 

4. The method of claim 1, wherein at least one of the calibration factors is a source 
location factor. 

5. The method of claim 1, wherein at least one of the calibration factors is a detector 
location factor. 



NY02:364961.1 



-43- 



A34927- 069225.0110 



6. The method of claim 1, wherein the medium has characteristics so as to highly 
scatter the radiation entering therein. 

7. The method of claim 1, wherein the radiation is an electromagnetic radiation. 

8. The method of claim 1, wherein the radiation is an infrared radiation. 

9. The method of claim 1, wherein the radiation comprises near-infrared photons. 

10. The method of claim 1 , wherein the radiation is a continuous-wave radiation. 

11. The method of claim 10, wherein at least one of the optical properties includes an 
absorption coefficient. 

12. The method of claim 10, wherein at least one of the optical properties includes a 
scattering coefficient. 

13. The method of claim 1, wherein the optical properties are spatially varying optical 
properties. 

14. The method of claim 1, wherein each of the optical properties is utilized to obtain 
amplitude and phase of the radiation exiting the medium, which are then used to 
determine the distribution. 
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15. The method of claim 1 , wherein each of the optical properties is utilized to obtain 
intensity and temporal delay of the radiation exiting the medium, which are then used to 
determine the distribution. 

16. The method of claim 1, further comprising the step of: 

d. before step (b), obtaining one or more parameters of the radiation received in 



17. The method of claim 16, wherein at least one of the parameters is fluence. 

18. The method of claim 17, wherein the step (b) comprises the substep of solving an 
inverse problem based on the one or more parameters and the one or more optical 
properties. 

19. The method of claim 18, wherein the solving substep comprises the substep of 
minimizing the equation: 



wherein index N max is a number of measurements of the received radiation, 
^Theory(x) is a theoretical fluence, 
$Meas(x) is a measured fluence, and 

x is a vector providing a value of a property distribution for each voxel 
representing a portion of the medium. 



step (a). 
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20. The method of claim 19, wherein the minimizing substep comprises utilizing a 
linear approximation. 

21 . The method of claim 20, wherein the linear approximation is a Rytov 
approximation. 

22. The method of claim 20, wherein the linear approximation is a Born 
approximation. 

23. The method of claim 16, wherein the obtained parameters are functions of the 
calibration factors. 

24. The method of claim 23, wherein at least one of the calibration factors is a source 
coupling factor. 

25. The method of claim 23, wherein at least one of the calibration factors is a 
detector coupling factor. 

26. The method of claim 23, wherein at least one of the calibration factors is a source 
location factor. 

27. The method of claim 23, wherein at least one of the calibration factors is a 
detector location factor. 

28. The method of claim 1, wherein at least one of the one or more optical properties 
includes an absorption coefficient value. 
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29. The method of claim 1 , wherein at least one of the one or more optical properties 
includes a scattering coefficient value. 

30. The method of claim 1, wherein the one or more optical properties include an 
absorption coefficient and a scattering coefficient. 

3 1 . The method of claim 1 , further comprising the step of: 

e. reconstructing at least one image of the medium by using the distribution 
determined in step (c). 

32. The method of claim 3 1 , further comprising the step of: 

f. displaying the reconstructed images for review by human operators. 

33. A method of measuring a distribution of at least one property within a medium, 
the method comprising: 

a. illuminating the medium with radiation from a plurality of sources; 

b. receiving radiation from the medium with a plurality of detectors; 

c. measuring at least one parameter of the received radiation; and 

d. calculating the distribution of the at least one property based on the at least 
one measured parameter by including source and detector calibration factors 
as freely varying quantities that are reconstructed in a model for a radiative 
transport within the medium. 
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34. The method of claim 33, wherein a probability that photons entering the medium 
will scatter greatly exceeds a probability that the photons entering the medium will be 
absorbed. 

35. The method of claim 33, wherein the radiation is an electromagnetic radiation. 

36. The method of claim 33, wherein the radiation is a continuous-wave radiation. 

37. The method of claim 33, wherein the radiation is an infrared radiation. 

38. The method of claim 33, wherein each source is spatially separated from each 
detector. 

39. The method of claim 33, wherein the measured parameters are at least one 
absorption parameter and at least one scattering parameter. 

40. The method of claim 33, wherein the measured parameters are at least one 
amplitude parameter and at least one phase parameter. 

41 . The method of claim 33, wherein the measured parameters are at least one 
amplitude parameter and at least one temporal off-set parameter. 

42. The method of claim 33, wherein the property distribution is a distribution of at 
least one of absorption coefficient values and scattering coefficient values. 
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43. The method of claim 33, further comprising the step of 

e. displaying the property distribution in at least one image. 

44. The method of claim 33, wherein the model includes a non-linear dependence on 
at least one measured parameter. 

45. The method of claim 33, wherein the model includes a linear dependence on at 
least one measured parameter. 

46. The method of claim 33, wherein the model includes a Rytov approximation. 

47. The method of claim 46, wherein the model includes the substep of minimizing 
the following equation: 



N max is a number of measurements of the received radiation, 
$Theory(x) is a theoretical fluence, 
^Meas(x) is a measured fluence, and 

x is a vector providing a value of the property distribution for each voxel of at 
least one image. 

48. The method of claim 33, wherein the model includes an arrangement for scaling 
all measured parameters taken as input to make the measured parameters dimensionless 
and to be of the same order as the source and detector calibration factors. 




wherein 
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49. A computer readable medium, comprising: 

a program, which when executed, is capable of causing a processor to: 

i. receive radiation exiting the medium, 

ii. derive one or more optical properties of the medium by using the received 
radiation and one or calibration factors, wherein the calibration factors are 
variables, and 

iii. determine as output the distribution of one or more optical properties 
based on the received radiation and one or more calibration factors. 

50. The computer readable medium of claim 49, wherein the program, when 
executed, is capable of causing the processor to reconstruct at least one image of the 
output. 

5 1 . The computer readable medium of claim 49, wherein at least one of the 
calibration factors is a source coupling factor. 

52. The computer readable medium of claim 49, wherein at least one of the 
calibration factors is a detector coupling factor. 

53. The computer readable medium of claim 49, wherein at least one of the 
calibration factors is a source location factor. 

54. The computer readable medium of claim 49, wherein at least one of the 
calibration factors is a detector location factor. 
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55. The computer readable medium of claim 49, wherein the medium has 
characteristics so as to highly scatter light entering therein. 

56. The computer readable medium of claim 49, wherein the optical properties are 
spatially varying optical properties. 

57. The computer readable medium of claim 49, wherein the program causes the 
processor to obtain one or more parameters of the radiation exiting the medium. 

58. The computer readable medium of claim 57, wherein at least one of the one or 
more parameters is fluence. 

59. The computer readable medium of claim 58, wherein the program, when 
executed, is capable of causing the processor to solve an inverse problem based on the 
one or more parameters and the one or more optical properties. 

60. The computer readable medium of claim 59, wherein the program, when 
executed, is capable of causing the processor to minimize the equation: 




wherein 



N max is a number of measurements of the received radiation, 



$Theory(x) is a theoretical fluence, 



$Meas(x) is a measured fluence, and 
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x is a vector providing a value of the property distribution for each voxel of at 
least one image. 

61 . The computer readable medium of claim 60, wherein the program, when 
executed, is capable of causing the processor to utilize a linear approximation for 
minimizing the equation F(x). 

62. The computer readable medium of claim 61 , wherein the linear approximation is a 
Born approximation. 

63. The computer readable medium of claim 61, wherein the linear approximation is a 
Rytov approximation. 

64. The computer readable medium of claim 63, wherein the obtained parameters of 
radiation exiting the medium are functions of the calibration factors. 

65 . A computer readable medium comprising: 

a program, which, when executed, is capable of causing a processor to: 



i. obtain at least one measured parameter of radiation exiting a medium, 

ii. calculate a distribution of at least one property based on the at least one 
measured parameter by including calibration factors as freely varying 
quantities that are reconstructed in a model for radiative transport within 
the medium, and 

iii. provide the property distribution as output. 
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66. The system of claim 65, wherein at least one of the calibration factors is a source 
coupling factor. 

67. The system of claim 65, wherein at least one of the calibration factors is a detector 
coupling factor. 

68. The system of claim 65, wherein at least one of the calibration factors is a source 
location factor. 

69. The system of claim 65, wherein at least one of the calibration factors is a detector 
location factor. 

70. The computer readable medium of claim 65, wherein the model includes a non- 
linear dependence on at least one measured parameter. 

71. The computer readable medium of claim 65, wherein the model includes a linear 
dependence on at least one measured parameter. 

72. The computer readable medium of claim 65, wherein the model includes a Rytov 
approximation. 
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73. The computer readable medium of claim 65, wherein the program, when 
executed, is capable of causing the processor to minimize the equation: 



wherein; 

N max is a number of measurements of the received radiation, 
^Theory(x) is a theoretical fluence, 
^Meas(x) is a measured fluence, and 

x is a vector providing a value of the property distribution for each voxel of at 
least one image. 

74. The computer readable medium of claim 73, wherein the model includes scaling 
all measured parameters taken as input to make the measured parameters dimensionless 
and of the same order as the source and detector calibration factors. 

75. A system for determining a distribution of one or more optical properties of a 
medium illuminated with radiation from one or more sources, comprising: 

a. one or more detectors for receiving radiation; and 

b. a processor coupled to the one or more detectors for deriving one or more 



optical properties from the radiation and one or more calibration factors, 
wherein the calibration factors are variables, and wherein the processor is 
programmed to determine the distribution from the derived optical properties. 
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76. The system of claim 75, further comprising: 



c. 



a computer database, operationally coupled to the processor, for storing a 



program that causes the processor to derive the optical properties and 



determine the distribution. 



77. The system of claim 76, wherein at least one of the calibration factors is a source 
coupling factor. 

78. The system of claim 76, wherein at least one of the calibration factors is a detector 
coupling factor. 

79. The system of claim 76, wherein at least one of the calibration factors is a source 
location factor. 

80. The system of claim 76, wherein at least one of the calibration factors is a detector 
location factor. 

81. The system of claim 76, wherein the medium has characteristics so as to highly 
scatter radiation entering therein. 

82. The system of claim 76, wherein the optical properties are spatially varying 
optical properties. 

83. The system of claim 76, wherein the radiation is an electromagnetic radiation. 

84. The system of claim 76, wherein the radiation is an infrared radiation. 
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85. The system of claim 76, wherein the radiation comprises near-infrared photons. 

86. The system of claim 76, wherein the radiation is a continuous-wave radiation. 

87. The system of claim 86, wherein at least one of the optical properties includes an 
absorption coefficient. 

88. The system of claim 86, wherein at least one of the optical properties includes a 
scattering coefficient. 

89. The system of claim 76, wherein each of the optical properties is utilized to obtain 
amplitude and phase of the radiation exiting the medium, which are then used to 
determine the distribution. 

90. The system of claim 76, wherein each of the optical properties is utilized to obtain 
intensity and temporal delay of the radiation exiting the medium, which are then used to 
determine the distribution. 

91 . The system of claim 76, wherein the one or more detectors are adapted to obtain 
one or more parameters of the radiation. 

92. The system of claim 91, wherein at least one of the parameters is fluence. 

93. The system of claim 92, wherein the processor is programmed to solve an inverse 
problem based on the one or more parameters and the one or more optical properties. 
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94. The system of claim 93, wherein the processor is programmed to minimize the 
equation: 



wherein: 

N max is a number of measurements of the received radiation, 
®Theory(x) is a theoretical fluence, 
^Meas(x) is a measured fluence, and 

x is a vector providing a value of the property distribution for each voxel of at 
least one image. 

95. The system of claim 94, wherein the processor is programmed to utilize a linear 
approximation. 

96. The system of claim 95, wherein the linear approximation is a Born 
approximation. 

97. The system of claim 95, wherein the linear approximation is a Rytov 
approximation. 

98. The system of claim 97, wherein the obtained parameters are functions of the 
calibration factors. 
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99. The system of claim 98, wherein at least one of the calibration factors is a source 
coupling factor. 

100. The system of claim 98, wherein at least one of the calibration factors is a detector 
coupling factor. 

101. The system of claim 98, wherein at least one of the calibration factors is a source 
location factor. 

102. The system of claim 98, wherein at least one of the calibration factors is a detector 
location factor. 

103. The system of claim 76, wherein at least one of the one or more optical properties 
includes an absorption coefficient. 

104. The system of claim 76, wherein at least one of the one or more optical properties 
includes a scattering coefficient. 

1 05. The system of claim 76, wherein the one or more optical properties include both 
absorption and scattering coefficients. 

106. The system of claim 76, wherein the processor is programmed to reconstruct at 
least one image of the medium by using the determined distribution. 

107. The system of claim 106, wherein the processor is programmed to display the 
reconstructed images for review by human operators. 
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108. A system for measuring a property distribution within a medium, the system 
comprising: 

a. a plurality of radiation sources; 

b. a plurality of detectors that convert the radiation received into signals; and 

c. a processor that processes the signals to provide values of the distribution of at 
least one property of the medium, wherein the values are obtained by applying 
a model for radiative transport in the medium that reconstructs calibration 
factors using freely varying quantities. 

109. The system of claim 108, wherein at least one of the calibration factors is a source 
coupling factor. 

110. The system of claim 108, wherein at least one of the calibration factors is a 
detector coupling factor. 

111. The system of claim 108, wherein at least one of the calibration factors is a source 
location factor. 

1 12. The system of claim 108, wherein at least one of the calibration factors is a 
detector location factor. 

113. The system of claim 108, wherein each source comprises an optical fiber and a 
laser, and wherein each detector comprises an optical fiber and a photodetector. 

114. The system of claim 108, wherein the radiation is an electromagnetic radiation. 



NY02:364961.1 



-59- 



• 



A34927- 069225.01 10 



115. The system of claim 108, wherein the radiation is a continuous-wave radiation. 

116. The system of claim 108, wherein the radiation is an infrared radiation. 

117. The system of claim 108, wherein each source is spatially separated from each 
detector. 

118. The system of claim 117, wherein each source extends in a first plane and each 
detector extends in a second plane which is opposite to the first plane. 

119. The system of claim 118, wherein each source is spatially separated from the 
other sources by about 2 centimeters and each detector is spatially separated from the 
other detectors by about 2 centimeters. 

120. The system of claim 108, wherein the processor provides output of at least one of 
absorption values and scattering values. 

121. The system of claim 108, further comprising a display to provide an image of at 
least one property of the medium. 

122. The system of claim 121, wherein the display is a computer screen. 

123. The system of claim 108, wherein the model is non-linear. 

124. The system of claim 108, wherein the model is linear. 

125. The system of claim 108, wherein the model includes a Rytov approximation. 
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126. The system of claim 125, wherein the model includes minimizing 

'"measr "l 

F(x) = J] [in (D neoryJ (x) - In O MeasJ J , 

i=l 

wherein: 

N max is a number of measurements of the received radiation, 
^Theory(x) is a theoretical fluence, 
$Meas(x) is a measured fluence, and 

x is a vector providing a value of the property distribution for each voxel of at 
least one image. 

127. The system of claim 126, wherein the model includes scaling parameters 
extracted from the signal as input to make them dimensionless and of the same order as 
the source and detector calibration factors. 

128. A software system which, when executed on a processing device, determines a 
distribution of one or more optical properties in a medium illuminated with radiation 
from one or more sources, the software system comprising: 

a processing subsystem which, when executed on the processing device, 
configures the processing device to perform the following: 

i. obtains parameters of the radiation exiting the medium, 

ii. derives one or more optical properties of the medium using one or more 
calibration factors and the obtained parameters, wherein the calibration 
factors are variables, and 
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iii. determines the distribution in the medium using the one or more derived 
optical properties. 

129. The software system of claim 128, wherein at least one of the calibration factors 
is a source coupling factor. 

130. The software system of claim 128, wherein at least one of the calibration factors 
is a detector coupling factor. 

131. The software system of claim 128, wherein at least one of the calibration factors 
is a source location factor. 

132. The software system of claim 128, wherein at least one of the calibration factors 
is a detector location factor. 

133. A system for determining a distribution of one or more optical properties of a 
medium illuminated with radiation from one or more sources, comprising: 

a. means for receiving radiation exiting the medium; 

b. means for deriving one or more optical properties of the medium using one or 
more calibration factors and the received radiation, wherein the calibration 
factors are variables; and 

c. means for determining the distribution in the medium using the derived 
optical properties. 
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134. The system of claim 133, wherein at least one of the calibration factors is a source 
coupling factor. 

135. The system of claim 133, wherein at least one of the calibration factors is a 
detector coupling factor. 

136. The system of claim 133, wherein at least one of the calibration factors is a source 
location factor. 

137. The system of claim 133, wherein at least one of the calibration factors is a 
detector location factor. 
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